Are WNL stars tracers of high metaUicity? 
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We present new atmosphere models for Wolf-Rayet stars that include a self-consistent solution 
of the wind hydrodynamics. We demonstrate that the formation of optically thick WR winds 
can be explained by radiative driving on Fe line opacities, implying a strong dependence on 
metaUicity {Z). Z-dependent model calculations for late-type WN stars show that these objects 
are very massive stars close to the Eddington limit, and that their formation is strongly favored 
for high metaUicity environments. 



1. PoWR hydrodynamic model atmospheres 

The Potsdam Wolf-Rayet (PoWR) hydrodynamic model atmospheres combine fully 
line-blanketed non-LTE models with the equations of hydrodynamics (for details see 
i Grafener fc Hamannl2005llHamann fc Grafeneil2003HKoesterke et" .l2002Hrxrafener et alJ 
l2002|) . The wind structure {p{r) and v{r)) and the temperature structure T(r) are com- 
puted in line with the full set of non-LTE populations, and the radiation field in the 
CO- moving frame (CMF). In contrast to all previous approaches, the radiative wind ac- 
celeration aiad is obtained by direct integration 

arad = - / Xi^^i^dz/, (1.1) 



instead of making use of the Sobolev approximation. In this way, complex processes 
like strong line overlap, or the redistribution of radiation, are automatically taken into 
account. Moreover, the models include small-scale wind clu mping (throughout this wor k 
we assume a clumping factor of Z? = 10, for details see also lHamann fc: Koesterkelll998|) . 
The models describe the conditions in WR atmospheres in a realistic manner, and provide 
synthetic spectra, i.e. they allow for a direct comparison with observations. 

Utilizing these models, we recently obtained the first fully self-c onsistent Wolf-Rayet 
wind model, for the case of an early-type WC star with strong lines ijGrafener fc HamannI 
Moreover, we hav e examined the mass loss from late-type WN stars and its de- 
pendence on metaUicity ijGrafener fc Hamannll2006at IGrafener fc Hamannll2006b|) . 



2. Spectral analyses of galactic WR stars 

A comprehensive study of galactic WR star s, based on sp ectral analyses with line- 
blanketed PoWR models llHamann et alJ l|2006l) for WN stars. iBarniske et all l|2006() for 
WC stars), revealed a bimodal WR subtype distribution in the HRD, where the H-rich 
WNL stars are located to the right of the ZAMS with luminosities above 10^ Lq, whereas 
the (mostly H-free) early- to intermediate WN subtypes, as well as the WC stars, show 
lower luminosities and hotter temperatures (see Fig.^. 

This dichotomy already implies that the H-rich WNL stars are the descendants of 
very massive stars, possibly still in the pase of central H-burning, whereas the earlier 
subtypes (including the WC stars) are more evolved, less massive. He-burning objects. 
Note, however, that distance estimates are only available for a small part of the WNL 
sample. Some of these objects thus might have lower luminosities and be the direct 
progenitors of the earlier subtypes. 
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Figure 1 . Rec e nt spe ctral analyses of R alactic WR stars with line-blanketed models, according to 
iHamann et alJ (|200fil ^ and lBarniske et al. (2006): symbols in light grey denote H-rich WR stars, 
whereas H-free objects are indicated in dark grey. WC stars are indicate d in black. For objects 
with large symbols distance estimates are available (|van der Huchtll200ir i. whereas objects with 
small symbols are calibrated by their spectral subtype. Evolutionary tracks for non-rotating 
massive stars (|Mevnct fc. Macdcr 20Q3) arc shown for comparison. 



3. Hydrodynamic atmosphere models for WNL stars 

In a recent work we have investigated the prop erties of th e luminous, H-rich WNL 
stars with our hydrodynamic PoWR models (.Grafener fc H amann 2006a|). The most 
important conclusion from that work is that WR-type mass loss is primarily triggered 
by high L/M ratios or, equivalently, Eddington factors Fe = Xoi*/47rcG'Mj, approaching 
unity. Note that high L /M ratios are expected for very massive stars and for He-burning 
objects, giving a natural explanation for the occurrence of the WR phenomenon. 

In Fig. 13 we show the results from grid computations for WNL stars with a fixed 
luminosity of 10^'^ Lq, and stellar temperatures in the range of 30-60 kK. For the 
stellar masses, values of 67 and 55 Mq are adopted, corresponding to Eddington fac- 
tors of Fc = 0.55 and 0.67. Notably, the mass loss strongly depends on Fc and T^,. The 
obtained synthetic spectra nicely reflect the observed sequence of weak-lined WNL sub- 
types, starting with WN 6 at 55 kK to WN 9 at 31 kK. From a more detailed investigation 
of the WN7 component in WR22, an eclipsing WR-l- binary system in Car OBI, we 
infer a stellar mass of 78 Mq (Fe — 0.67), in agreement iRauw et al. (199^ who obtained 
72±3 from the binary orbit. Such high stellar masses imply that the weak-lined WNL 
stars are still in the pase of central H-burning, suggesting an evolutionary sequence of 
the form O WNL LBV WN WC for very massive stars. 
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Figure 2. Wind models for galactic WNL stars: mass loss rates for different stellar temperatures 
and Eddington factors To. The corresponding spectral subtypes are indicated together with 
WR numbers of specific galactic objects (according to van dcr Hucht 200^, in brackets), showing 
a good agreement with the synthetic line spectra. Note that the models are computed for a fixed 
stellar luminosity of lO**'^ Lq. The WN 7 model (WR22) is slightly offset from the standard grid 
models because it is calculated with an enhanced hydrogen abundance (see text). 

4. WNL stars at different metallicities 

In Fig.Elw e present results from a grid of models for luminous WNL stars at different 
metallicities (iGrafener fc Hamannll2006a|) . The models are computed for a fixed lumi- 
nosity of 10^ '^ Lq, a stellar temperature of Tj, = 45 kK, and a hydrogen surface mass 
fraction of = 0.4. In addition to the metal abundances we have varied the Eddington 
factor Fe (or equivalently the stellar mass). Note that we have scaled all metals with Z, 
assuming a CNO-processed (i.e. N-enriched) WN surface composition. The wind driving 
in our mod els is chiefly due t o radi ation pressure on Fe-group line opacities. In accor- 
dance with IVink fc de KoteJ l)2005|) we thus find a strong dependence on Z. However, 
as in our previous computations, the proximity to the Eddington limit plays an equally 
important role. We find that optically thick winds with high WR-type mass loss rates 
are formed over the whole range of metallicities, from 1/1000-3 Zq, if the stars get close 
enough to the Eddington limit. Only the limiting value of Fg where the WR-type winds 
start to form, changes. For solar Z, values of Fg = 0.5-0.6 are leading to the formation 
of weak-lined WNL stars. As we have seen for the case of WR22, this corresponds to 
very massive, slightly over-luminous stars in a late phase of H-burning. Note that it is 
indeed observed that the most massive stars in very you ng galactic clusters ar e in the 
WNL phase fe.ef..'Figer et al.V2002').'N aiarro et alJ ^|2004^ for the Arches cluster: IPrissenl 
1.1999), Crowther & Dcssart (1998) for NGC3603). 

For higher values of Z, the limit for the formation of WR-type winds shifts towards 
even lower values of Fg. For 3 Zq, we find that stars with Fe Pi 0.4 already show typical 
WNL mass loss rates (see Fig.OJ. This corresponds to objects with 120 Mq on the ZAMS. 
We thus expect that metal rich stars with very high masses already start their life in 
the WNL phase, i.e., the occurrence of WNL stars in young massive clusters is strongly 
favored for high metallicities. 
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Figure 3. WNL star mass loss over a broad range of metallicities (Z): Mass loss rates (top) 
and terminal wind velocities (bottom), as obtained from our hydrodynamic models, are plotted 
vs. the Eddington factor Fg. The soUd curves indicate model series for WNL stars where Fe is 
varied for a given value of Z. The dashed grey lines indicate models with 3Zq and Xu = 0.7, 
corresponding to very massive, metal rich stars on the ZAMS. 
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